afs_orig = read.csv("data/MGRB_GnomAD_SweGen_cancer_UKBB_AFs_outerjoined_ss-auto-gwas-snps_hcr.csv.xz", stringsAsFactors = FALSE, header = TRUE)
models_orig = read.csv("data/manual_polygenic_scores.hcr_tag_rescued.csv", stringsAsFactors = FALSE, header = TRUE)
# Key allele frequencies by VID
afs_orig$vid = paste(afs_orig$chrom, afs_orig$pos, afs_orig$ref, afs_orig$alt, sep = ":")
afs_orig = afs_orig[,!(colnames(afs_orig) %in% c("chrom", "pos", "ref", "alt"))]
# Add gnomad AFs to the models for imputation of missing variants
temp.gnomad_af = (afs_orig$nAA_gnomad*2 + afs_orig$nRA_gnomad) / (2*(afs_orig$nAA_gnomad + afs_orig$nRA_gnomad + afs_orig$nRR_gnomad))
models_orig$aaf = temp.gnomad_af[match(models_orig$vid, afs_orig$vid)]
# Create a full UKBB cohort by combining the age-stratified numbers
afs_orig$nRR_ukbb = afs_orig$nRR_ukbb_0_55 + afs_orig$nRR_ukbb_55_60 + afs_orig$nRR_ukbb_60_65 + afs_orig$nRR_ukbb_65_70 + afs_orig$nRR_ukbb_70_75 + afs_orig$nRR_ukbb_75_inf
afs_orig$nRA_ukbb = afs_orig$nRA_ukbb_0_55 + afs_orig$nRA_ukbb_55_60 + afs_orig$nRA_ukbb_60_65 + afs_orig$nRA_ukbb_65_70 + afs_orig$nRA_ukbb_70_75 + afs_orig$nRA_ukbb_75_inf
afs_orig$nAA_ukbb = afs_orig$nAA_ukbb_0_55 + afs_orig$nAA_ukbb_55_60 + afs_orig$nAA_ukbb_60_65 + afs_orig$nAA_ukbb_65_70 + afs_orig$nAA_ukbb_70_75 + afs_orig$nAA_ukbb_75_inf
afs_orig$nmissing_ukbb = afs_orig$nmissing_ukbb_0_55 + afs_orig$nmissing_ukbb_55_60 + afs_orig$nmissing_ukbb_60_65 + afs_orig$nmissing_ukbb_65_70 + afs_orig$nmissing_ukbb_70_75 + afs_orig$nmissing_ukbb_75_inf
# Convert afs from wide to long format
library(reshape2)
afs_long = melt(afs_orig, id.vars = c("rsid", "negstrand", "vid"), value.name = "count")
afs_long$cohort = gsub("^n(RR|RA|AA|missing)_", "", afs_long$variable)
afs_long$variable = gsub("_.*", "", afs_long$variable)
afs = dcast(afs_long, vid + rsid + negstrand + cohort ~ variable, value.var = "count")
afs = afs[,c("vid", "rsid", "negstrand", "cohort", "nRR", "nRA", "nAA", "nmissing")]
# Exclude ASRB samples -- prelim examination suggests they are rather
# poor quality, and we are not interested in their PRS distributions
# anyway. Also exclude the various MGRB filtration options, as they
# apply only to rare variants. Exclude SweGen as we don't have a good
# HQ bed for it.
cohorts.sel = c(
"mgrborig",
"gnomad",
"ukbb", "ukbb_0_55", "ukbb_55_60", "ukbb_60_65", "ukbb_65_70", "ukbb_70_75", "ukbb_75_inf",
"cancercrc", "cancermel", "cancernms", "cancerbrca_f", "cancerpca_m",
"nocancer", "nocancer_f", "nocancer_m", "mgrborig_f", "mgrborig_m")
cohorts.main = c("mgrborig", "gnomad", "ukbb")
afs = afs[afs$cohort %in% cohorts.sel,]
Choose polygenic models with at least 10 loci, with the exception of ShortLifespan:Deelen:10.1093/hmg/ddu139 (only six loci passing filters). For cancers, choose polygenic models only for cancers with a positive control cohort. In the case of multiple models for the same disorder, choose the most recent original publication where possible (ie exclude “meta” signatures if a good original report is available).
# Excluded due to insufficient size:
# "CancerOfBladder:Fritsche:10.1016/j.ajhg.2018.04.001",
# "LymphoidLeukemiaAcute:Fritsche:10.1016/j.ajhg.2018.04.001",
# "LymphoidLeukemiaChronic:Fritsche:10.1016/j.ajhg.2018.04.001",
# "MalignantNeoplasmOfTestis:Fritsche:10.1016/j.ajhg.2018.04.001",
# "NonHodgkinsLymphoma:Fritsche:10.1016/j.ajhg.2018.04.001",
# "PancreaticCancer:Fritsche:10.1016/j.ajhg.2018.04.001",
# "APOE_rs429358:NA:NA",
# Excluded because a better alternative was available
# "CancerOfProstate:Fritsche:10.1016/j.ajhg.2018.04.001",
# "ColorectalCancer:Fritsche:10.1016/j.ajhg.2018.04.001",
# "MelanomasOfSkin:Fritsche:10.1016/j.ajhg.2018.04.001",
# "BreastCancer:Li:10.1038/gim.2016.43",
# "BreastCancerFemale:Fritsche:10.1016/j.ajhg.2018.04.001",
# Excluded because of issues with population-specific alleles between UK and European popns
# "BasalCellCarcinoma:Chahal:10.1038/ncomms12510",
# "BasalCellCarcinoma:Fritsche:10.1016/j.ajhg.2018.04.001",
models.sel = c(
"AF:Lubitz:10.1161/CIRCULATIONAHA.116.024143",
"DiastolicBP:Warren:10.1038/ng.3768",
"EOCAD:Theriault:10.1161/circgen.117.001849",
"PulsePressure:Warren:10.1038/ng.3768",
"SystolicBP:Warren:10.1038/ng.3768",
"AlzheimersDisease:Lambert:10.1038/ng.2802",
"ShortLifespan:Deelen:10.1093/hmg/ddu139",
"Height:Wood:10.1038/ng.3097",
"BreastCancer:Michailidou:10.1038/nature24284",
"ColorectalCancer:Schumacher:10.1038/ncomms8138",
"Melanoma:Law:10.1038/ng.3373",
"ProstateCancer:Hoffmann:10.1158/2159-8290.CD-15-0315"
)
models = models_orig[models_orig$id %in% models.sel,]
# Drop variants with low genotyping rate in any cohort in which that
# variant was detected. Use the threshold of 97% genotyping rate
library(plyr)
temp.lowgt = ddply(afs[!is.na(afs$nRR),], .(cohort), function(d) mean(d$nmissing / (d$nRR + d$nRA + d$nAA + d$nmissing) >= 0.03))
temp.lowgt
## cohort V1
## 1 cancerbrca_f 1.637465e-04
## 2 cancercrc 1.228099e-04
## 3 cancermel 1.228099e-04
## 4 cancernms 8.187326e-05
## 5 cancerpca_m 4.093663e-05
## 6 gnomad 1.598146e-03
## 7 mgrborig 8.074609e-05
## 8 mgrborig_f 8.074609e-05
## 9 mgrborig_m 8.074609e-05
## 10 nocancer 8.187326e-05
## 11 nocancer_f 4.093663e-05
## 12 nocancer_m 1.228099e-04
## 13 ukbb 5.066030e-02
## 14 ukbb_0_55 5.120115e-02
## 15 ukbb_55_60 5.097580e-02
## 16 ukbb_60_65 5.097580e-02
## 17 ukbb_65_70 5.057015e-02
## 18 ukbb_70_75 5.038987e-02
## 19 ukbb_75_inf 4.989408e-02
# The UKBB samples have rather a lot of dropouts here: ~ 5% of loci
# have a gt rate under 97%.
temp.gt_rate = 1 - daply(afs[!is.na(afs$nRR),], .(vid), function(d) max(d$nmissing / (d$nRR + d$nRA + d$nAA + d$nmissing)))
mean(temp.gt_rate < 0.97)
## [1] 0.04893795
afs = afs[afs$vid %in% names(temp.gt_rate)[temp.gt_rate >= 0.97],]
mean(models$vid %in% names(temp.gt_rate)[temp.gt_rate >= 0.97])
## [1] 0.9484241
# ~5.2% of model loci lost by this filter
models = models[models$vid %in% names(temp.gt_rate)[temp.gt_rate >= 0.97],]
# Create a set of AFs for variants that have VCF entries in every cohort.
# Note that given the relatively small size of some cohorts, this tends to
# preferentially exclude rare variants from consideration, and will probably
# attenuate the GRS differences.
afs.nmissing_per_cohort = tapply(is.na(afs$nRR), afs$vid, sum)
afs.nomissing = afs[!(afs$vid %in% names(afs.nmissing_per_cohort[afs.nmissing_per_cohort > 0])),]
nrow(afs.nomissing) / nrow(afs)
## [1] 0.8663734
mean(models$vid %in% afs.nomissing$vid)
## [1] 0.6628399
Note that we examine only loci in the PRSes.
temp.model_vids = intersect(afs.nomissing$vid, models$vid)
temp.afs = afs.nomissing[afs.nomissing$vid %in% temp.model_vids,]
temp.afs$fA = (temp.afs$nAA*2 + temp.afs$nRA) / (2*(temp.afs$nAA + temp.afs$nRA + temp.afs$nRR))
temp.afs = acast(temp.afs, vid ~ cohort, value.var = "fA", fill = NA)
pairs(temp.afs[,cohorts.main], pch = ".")
temp.overall_cohort.pvals = sapply(cohorts.main[1:(length(cohorts.main)-1)], function(cohort_1) {
cohort_1_idx = which(cohorts.main == cohort_1)
afs_1 = temp.afs[,cohort_1]
sapply(cohorts.main[(cohort_1_idx+1):length(cohorts.main)], function(cohort_2) {
afs_2 = temp.afs[,cohort_2]
test = wilcox.test(afs_1 - afs_2)
test$p.value
})
})
temp.overall_cohort.pvals
## $mgrborig
## gnomad ukbb
## 0.69792760 0.01866417
##
## $gnomad
## ukbb
## 0.3782686
p.adjust(unlist(temp.overall_cohort.pvals), "holm")
## mgrborig.gnomad mgrborig.ukbb gnomad.ukbb
## 0.75653715 0.05599252 0.75653715
g.test = function(tbl)
{
expected = outer(rowSums(tbl), colSums(tbl)) / sum(tbl)
logoe = log(tbl/expected)
logoe[tbl == 0] = 0
stat = 2*sum(tbl*logoe)
pchisq(stat, prod(dim(tbl)-1), lower.tail = FALSE)
}
temp.locus_cohort.pvals = ddply(afs.nomissing[afs.nomissing$vid %in% temp.model_vids & afs.nomissing$cohort %in% cohorts.main,], .(vid), function(d) {
nR = d$nRR*2 + d$nRA
nA = d$nAA*2 + d$nRA
g.test(cbind(nR, nA))
# g.test(as.matrix(d[,c("nRR", "nRA", "nAA")]))
})
colnames(temp.locus_cohort.pvals)[2] = "p.value"
temp.locus_cohort.pvals$p.value.mtc = p.adjust(temp.locus_cohort.pvals$p.value, "BH")
# temp.locus_cohort.pvals$p.value.mtc is now calibrated for average false rejection rate
# (ie calling a SNP population-associated when it in fact isn't).
mean(temp.locus_cohort.pvals$p.value.mtc < 0.01)
## [1] 0.4721461
temp.locus_cohort.maxdeltaaf = ddply(afs.nomissing[afs.nomissing$vid %in% temp.model_vids & afs.nomissing$cohort %in% cohorts.main,], .(vid), function(d) {
nR = d$nRR*2 + d$nRA
nA = d$nAA*2 + d$nRA
AAF = nA/(nA+nR)
max(AAF) - min(AAF)
})
colnames(temp.locus_cohort.maxdeltaaf)[2] = "maxdeltaaaf"
temp.locus_cohort.maxdeltaaf[order(temp.locus_cohort.maxdeltaaf$maxdeltaaaf),]
## vid maxdeltaaaf
## 629 20:41851935:G:A 4.804901e-05
## 175 11:244552:A:G 5.472387e-04
## 680 3:163838015:A:C 5.854864e-04
## 832 5:176517326:T:C 6.520533e-04
## 849 5:58337481:T:G 6.896540e-04
## 576 2:232268312:T:C 8.648204e-04
## 64 1:242034263:A:G 8.926576e-04
## 557 2:19942473:G:A 8.944422e-04
## 783 4:61995613:A:G 9.804085e-04
## 343 15:32993111:C:T 1.013646e-03
## 125 10:12943973:C:T 1.036166e-03
## 428 17:43216281:C:T 1.120399e-03
## 923 6:43711981:T:C 1.242711e-03
## 87 1:65010606:T:G 1.410640e-03
## 1076 9:22062134:G:T 1.487932e-03
## 943 6:81792063:G:T 1.508864e-03
## 368 15:72842705:G:A 1.511067e-03
## 716 3:62133492:G:A 1.517305e-03
## 754 4:147993702:A:G 1.550150e-03
## 49 1:215046892:G:A 1.631525e-03
## 45 1:203766331:A:G 1.651584e-03
## 775 4:39503196:A:G 1.678110e-03
## 799 5:108625324:C:A 1.746323e-03
## 34 1:17308254:T:C 1.809376e-03
## 462 18:19450303:A:G 1.867013e-03
## 460 17:800593:T:C 1.957503e-03
## 661 21:45867411:G:A 1.968592e-03
## 290 13:55934157:A:G 1.984521e-03
## 1095 9:99203606:T:C 2.015219e-03
## 543 2:179786068:T:C 2.214285e-03
## 221 12:125441159:T:C 2.246113e-03
## 364 15:67455630:C:T 2.312045e-03
## 444 17:64783539:C:T 2.316572e-03
## 1064 9:118826916:G:A 2.361107e-03
## 383 15:98615560:C:T 2.370182e-03
## 698 3:190815978:A:G 2.376253e-03
## 17 1:151259043:C:T 2.419978e-03
## 950 7:107259721:T:C 2.424953e-03
## 918 6:36659932:C:T 2.450101e-03
## 670 3:136107549:G:A 2.505036e-03
## 13 1:14105298:G:A 2.546910e-03
## 871 6:118569679:T:G 2.552647e-03
## 993 8:110115372:C:T 2.598390e-03
## 158 11:11563879:C:T 2.615930e-03
## 300 13:80618435:A:G 2.616690e-03
## 965 7:148629759:C:T 2.621543e-03
## 701 3:27416013:C:T 2.622335e-03
## 548 2:191227755:A:G 2.638522e-03
## 477 18:57323149:C:A 2.650551e-03
## 1005 8:128106880:A:C 2.808417e-03
## 36 1:177279412:G:A 2.814433e-03
## 791 4:82204091:A:G 2.816117e-03
## 936 6:75452066:T:C 2.865091e-03
## 502 19:42683964:C:T 2.865981e-03
## 801 5:113748571:C:T 2.886389e-03
## 633 20:50141264:T:C 2.900362e-03
## 883 6:149608874:G:A 2.955662e-03
## 995 8:115698881:G:A 2.962404e-03
## 482 18:77222862:T:G 2.971526e-03
## 332 14:90636206:G:A 3.008568e-03
## 572 2:219903723:C:T 3.051369e-03
## 650 21:27208935:G:T 3.062235e-03
## 827 5:171285632:C:T 3.075935e-03
## 606 2:56113538:G:A 3.098396e-03
## 812 5:141681788:G:A 3.104595e-03
## 1067 9:127900996:T:C 3.347011e-03
## 687 3:178467852:A:G 3.415386e-03
## 618 20:19465907:G:A 3.418120e-03
## 691 3:185313855:A:G 3.418663e-03
## 153 10:97805074:G:A 3.442529e-03
## 942 6:81038921:T:G 3.447571e-03
## 994 8:110581794:A:G 3.470129e-03
## 180 11:45643843:C:T 3.483223e-03
## 466 18:29088958:C:T 3.498834e-03
## 865 6:109285189:A:G 3.553083e-03
## 390 16:2345388:T:C 3.579194e-03
## 58 1:227191011:A:G 3.623518e-03
## 777 4:42308930:C:T 3.649195e-03
## 600 2:44392271:C:A 3.688753e-03
## 939 6:7709052:G:T 3.701351e-03
## 723 3:89530057:C:A 3.732115e-03
## 739 4:111765495:C:T 3.739535e-03
## 1000 8:123980551:C:T 3.752414e-03
## 1038 8:57095808:T:C 3.768970e-03
## 808 5:134499092:C:A 3.775173e-03
## 790 4:82184049:C:T 3.790315e-03
## 342 15:30167418:G:A 3.798314e-03
## 8 1:11852516:C:T 3.841580e-03
## 674 3:141163045:T:C 3.905120e-03
## 686 3:178343683:C:T 3.911620e-03
## 362 15:66997087:A:G 3.913701e-03
## 96 1:89360520:T:G 3.941459e-03
## 1085 9:92498089:C:T 3.992649e-03
## 437 17:56431549:A:G 4.036653e-03
## 744 4:119121695:C:T 4.041940e-03
## 480 18:74983055:A:G 4.072636e-03
## 164 11:128478885:C:A 4.108551e-03
## 82 1:54119578:G:A 4.159979e-03
## 1037 8:56998480:G:A 4.187366e-03
## 1080 9:6365683:A:C 4.206827e-03
## 433 17:47944836:G:A 4.216069e-03
## 571 2:219154781:G:A 4.227556e-03
## 1066 9:126094366:T:C 4.254321e-03
## 438 17:59483766:C:T 4.254752e-03
## 647 21:16520832:G:A 4.264179e-03
## 212 12:114804580:A:C 4.266542e-03
## 268 13:111100780:C:T 4.279067e-03
## 578 2:232796610:T:C 4.294600e-03
## 315 14:53391680:A:G 4.313766e-03
## 307 14:100599519:T:C 4.321819e-03
## 703 3:38036914:C:T 4.336252e-03
## 848 5:56254485:G:T 4.379536e-03
## 919 6:39134099:T:C 4.392131e-03
## 699 3:191111160:A:G 4.457648e-03
## 73 1:32371442:C:T 4.488190e-03
## 901 6:17589375:T:C 4.528051e-03
## 40 1:195866450:T:G 4.531433e-03
## 440 17:62161933:A:G 4.543676e-03
## 176 11:2810731:C:T 4.564121e-03
## 327 14:70609793:G:A 4.566616e-03
## 309 14:24830850:T:C 4.616386e-03
## 273 13:22453087:T:C 4.668375e-03
## 467 18:29363000:G:A 4.700294e-03
## 337 14:93784292:A:C 4.705838e-03
## 930 6:51450906:C:T 4.742416e-03
## 114 10:114773927:A:G 4.743066e-03
## 831 5:175947118:G:T 4.774634e-03
## 539 2:177557020:A:G 4.794079e-03
## 299 13:78474468:A:C 4.812120e-03
## 43 1:20030836:G:A 4.816278e-03
## 174 11:23197362:G:A 4.825182e-03
## 83 1:56965664:C:T 4.866564e-03
## 510 19:56323209:T:C 4.882986e-03
## 395 16:52599188:C:T 4.899116e-03
## 135 10:44777188:C:A 4.900737e-03
## 839 5:32821168:A:G 4.908845e-03
## 1010 8:129194641:C:T 4.921311e-03
## 134 10:44480811:T:G 4.943897e-03
## 267 13:111049623:T:C 4.952131e-03
## 170 11:2022804:A:G 4.978664e-03
## 414 17:1353920:T:G 5.093687e-03
## 733 4:109408608:G:A 5.102530e-03
## 574 2:221397008:T:C 5.228419e-03
## 426 17:40565926:T:G 5.228988e-03
## 249 12:6385727:C:T 5.238769e-03
## 722 3:74522679:G:A 5.286244e-03
## 523 2:142077308:T:C 5.309499e-03
## 347 15:39531658:A:G 5.312263e-03
## 351 15:51269629:A:G 5.343431e-03
## 384 15:99219598:C:T 5.350541e-03
## 348 15:39597034:T:C 5.358179e-03
## 264 12:96027759:A:G 5.362742e-03
## 999 8:117630683:A:C 5.369172e-03
## 810 5:137044526:T:C 5.376612e-03
## 1059 9:114850255:A:C 5.388091e-03
## 770 4:29571133:C:T 5.413917e-03
## 376 15:81836638:G:A 5.443466e-03
## 90 1:77559860:G:A 5.447211e-03
## 755 4:148400819:C:A 5.449392e-03
## 768 4:184529029:T:G 5.475820e-03
## 132 10:35150975:A:G 5.499058e-03
## 1073 9:139111870:G:T 5.525798e-03
## 152 10:96632253:T:C 5.534003e-03
## 542 2:179349227:T:G 5.541236e-03
## 142 10:75683793:G:A 5.549552e-03
## 835 5:31540303:A:G 5.551595e-03
## 880 6:13722523:G:A 5.632665e-03
## 103 10:102604514:G:A 5.670602e-03
## 341 15:101762539:C:T 5.709396e-03
## 50 1:21584585:A:G 5.722582e-03
## 141 10:75416789:A:G 5.731772e-03
## 210 12:114556438:C:T 5.756619e-03
## 814 5:142928348:C:A 5.757560e-03
## 1069 9:133464084:A:G 5.818366e-03
## 676 3:14928077:G:T 5.820147e-03
## 945 6:82442022:G:A 5.842640e-03
## 429 17:46096276:C:T 5.845469e-03
## 547 2:190404135:C:T 5.887918e-03
## 457 17:7571752:T:G 5.938427e-03
## 764 4:175846426:C:A 5.946180e-03
## 913 6:34195011:G:T 5.966954e-03
## 561 2:207948768:C:T 5.969859e-03
## 261 12:93919840:C:T 6.005322e-03
## 350 15:48914926:A:G 6.033136e-03
## 653 21:32487917:C:T 6.037518e-03
## 77 1:41101572:A:G 6.077944e-03
## 521 2:134434824:T:C 6.094110e-03
## 208 12:104344836:A:G 6.112237e-03
## 644 20:6572014:A:C 6.117750e-03
## 663 3:111560228:G:A 6.119474e-03
## 329 14:74990746:A:G 6.133487e-03
## 333 14:91841069:A:G 6.134627e-03
## 756 4:148763687:C:T 6.140946e-03
## 788 4:7912333:G:A 6.149031e-03
## 41 1:196059593:T:C 6.159112e-03
## 494 19:19591066:A:G 6.165739e-03
## 785 4:71267096:G:A 6.208946e-03
## 636 20:55015166:A:C 6.210536e-03
## 992 8:101677386:T:C 6.243719e-03
## 735 4:111662573:G:T 6.245742e-03
## 105 10:104341435:T:C 6.264502e-03
## 416 17:29234505:C:T 6.265179e-03
## 616 20:17771113:A:G 6.270317e-03
## 1035 8:53978574:C:T 6.281451e-03
## 245 12:57527283:T:C 6.284469e-03
## 196 11:69163161:C:T 6.351210e-03
## 1060 9:117011595:C:T 6.355275e-03
## 690 3:184513018:T:C 6.371485e-03
## 119 10:123093901:C:T 6.395301e-03
## 580 2:241882785:C:T 6.429689e-03
## 677 3:154711854:C:T 6.459046e-03
## 99 1:95297449:G:A 6.465789e-03
## 475 18:53352082:T:G 6.471711e-03
## 349 15:41790987:C:T 6.483439e-03
## 968 7:16984280:C:T 6.522927e-03
## 607 2:564079:G:T 6.541430e-03
## 527 2:145258445:C:T 6.554976e-03
## 286 13:50718468:C:T 6.593988e-03
## 107 10:10499796:T:G 6.613064e-03
## 982 7:40447971:G:A 6.623932e-03
## 798 5:108113344:C:T 6.648022e-03
## 1054 9:109518208:A:G 6.686976e-03
## 311 14:28338474:G:T 6.691137e-03
## 276 13:27665405:C:T 6.691232e-03
## 171 11:21448255:G:A 6.708041e-03
## 555 2:198711900:C:A 6.710727e-03
## 750 4:138229547:C:T 6.716808e-03
## 575 2:225030129:T:C 6.726515e-03
## 328 14:73932966:T:C 6.754591e-03
## 765 4:17936634:T:C 6.756833e-03
## 289 13:55683505:G:T 6.762488e-03
## 797 5:108073934:C:A 6.763213e-03
## 511 19:7184762:A:G 6.764992e-03
## 774 4:39343940:C:A 6.784044e-03
## 63 1:237146379:A:G 6.848351e-03
## 935 6:75428815:C:T 6.877728e-03
## 567 2:218146080:C:T 6.879259e-03
## 182 11:56230069:T:C 6.883608e-03
## 675 3:142535505:C:T 6.885987e-03
## 503 19:44082429:C:A 6.898736e-03
## 88 1:67510474:A:G 6.916826e-03
## 464 18:2766938:A:G 6.954013e-03
## 984 7:55889334:G:A 6.954748e-03
## 899 6:166329862:C:T 6.963918e-03
## 833 5:176675423:A:G 6.965767e-03
## 875 6:122528771:C:T 6.968636e-03
## 837 5:32694942:C:T 6.995458e-03
## 1055 9:109599046:G:A 6.999378e-03
## 214 12:115836522:A:G 7.002404e-03
## 541 2:178685534:A:G 7.016063e-03
## 941 6:80956208:C:T 7.040035e-03
## 811 5:137419989:C:T 7.057785e-03
## 417 17:29247715:G:A 7.081359e-03
## 667 3:129050943:A:G 7.091920e-03
## 591 2:37438764:A:G 7.093670e-03
## 324 14:68660428:T:C 7.147187e-03
## 288 13:51106788:G:A 7.163935e-03
## 425 17:38632852:T:C 7.193085e-03
## 544 2:179815641:A:G 7.219288e-03
## 137 10:53532468:T:C 7.219707e-03
## 1056 9:110837073:A:G 7.223326e-03
## 30 1:169591164:T:C 7.288344e-03
## 94 1:86330770:C:T 7.303259e-03
## 681 3:16670272:C:T 7.324530e-03
## 967 7:150594014:C:T 7.326090e-03
## 151 10:91002927:C:T 7.341830e-03
## 778 4:4308513:T:C 7.358291e-03
## 515 2:10566804:A:G 7.364454e-03
## 400 16:72079657:C:T 7.395560e-03
## 1026 8:25298895:A:G 7.402153e-03
## 802 5:11638501:A:G 7.433651e-03
## 669 3:13555836:C:T 7.456637e-03
## 817 5:15447939:C:T 7.460525e-03
## 763 4:1693931:C:T 7.496437e-03
## 781 4:56399648:T:C 7.519776e-03
## 281 13:33723244:C:T 7.547754e-03
## 1081 9:78542286:G:A 7.550492e-03
## 1089 9:96893945:A:G 7.579373e-03
## 470 18:45888770:G:A 7.594407e-03
## 79 1:46024454:G:A 7.615989e-03
## 895 6:161108536:C:T 7.630865e-03
## 101 10:101345366:C:T 7.642579e-03
## 394 16:51094038:C:T 7.672765e-03
## 154 11:103673277:A:C 7.677657e-03
## 254 12:67872464:C:T 7.683610e-03
## 257 12:78839926:T:C 7.688332e-03
## 200 11:74345550:T:G 7.706280e-03
## 352 15:51530495:G:A 7.710161e-03
## 886 6:152157881:C:T 7.715515e-03
## 902 6:18240002:T:C 7.723580e-03
## 500 19:39592303:C:T 7.727099e-03
## 657 21:36481863:C:A 7.741424e-03
## 338 14:98587630:T:C 7.745059e-03
## 569 2:218616633:C:T 7.746794e-03
## 468 18:32395456:G:A 7.802055e-03
## 413 17:12449192:A:G 7.803015e-03
## 590 2:36768875:G:A 7.830523e-03
## 864 6:105392745:C:T 7.837228e-03
## 388 16:20880422:C:T 7.837748e-03
## 531 2:163284067:T:C 7.838424e-03
## 260 12:90252500:G:A 7.866925e-03
## 84 1:57015668:A:G 7.873771e-03
## 181 11:47557871:T:C 7.885528e-03
## 197 11:69367118:A:C 7.920853e-03
## 21 1:154714006:G:A 7.932151e-03
## 38 1:184007119:C:T 7.938685e-03
## 439 17:61936147:C:T 7.946537e-03
## 248 12:59956923:A:G 7.976395e-03
## 598 2:43541386:G:A 8.143734e-03
## 929 6:50615935:G:A 8.173286e-03
## 301 13:80705315:G:T 8.197344e-03
## 558 2:202181247:C:T 8.198081e-03
## 244 12:56680636:A:G 8.199471e-03
## 186 11:61571478:T:C 8.214525e-03
## 779 4:48687351:A:G 8.217647e-03
## 593 2:37765176:A:G 8.225633e-03
## 1002 8:128011937:T:C 8.234045e-03
## 759 4:156513518:T:C 8.264346e-03
## 1084 9:89108161:C:T 8.265327e-03
## 242 12:53273904:G:A 8.273572e-03
## 922 6:41924931:A:G 8.282893e-03
## 61 1:230188748:G:A 8.332156e-03
## 1022 8:17338553:C:T 8.340575e-03
## 147 10:80942721:A:G 8.346950e-03
## 308 14:24764114:A:G 8.378492e-03
## 627 20:40964651:C:T 8.383925e-03
## 666 3:11647173:A:G 8.395610e-03
## 780 4:5016883:G:A 8.427405e-03
## 98 1:9441949:C:A 8.430541e-03
## 822 5:166392055:C:T 8.458541e-03
## 334 14:92427348:T:C 8.459718e-03
## 236 12:40952895:G:A 8.470985e-03
## 95 1:87633477:C:T 8.472516e-03
## 411 16:990815:C:T 8.497403e-03
## 975 7:25871109:C:T 8.499375e-03
## 536 2:173311553:A:G 8.515985e-03
## 1003 8:128093297:C:T 8.536973e-03
## 917 6:36622900:C:A 8.563625e-03
## 160 11:120257495:A:G 8.589033e-03
## 602 2:46921285:G:A 8.599505e-03
## 44 1:201872264:C:T 8.603487e-03
## 15 1:146692373:C:T 8.633313e-03
## 292 13:60917626:C:T 8.637416e-03
## 720 3:67416322:T:C 8.638049e-03
## 187 11:61597972:A:G 8.644408e-03
## 487 19:12186611:A:G 8.650756e-03
## 216 12:117365506:G:A 8.710141e-03
## 845 5:44706498:A:G 8.749387e-03
## 253 12:67089330:G:A 8.753526e-03
## 911 6:32855458:T:C 8.763688e-03
## 638 20:59578269:T:C 8.772960e-03
## 525 2:144159905:A:C 8.783298e-03
## 456 17:74574826:C:T 8.814369e-03
## 1058 9:113792706:C:A 8.840161e-03
## 858 5:80391253:A:G 8.857158e-03
## 1065 9:124412948:T:G 8.863803e-03
## 220 12:123822711:G:A 8.906989e-03
## 1082 9:79960141:G:A 8.913657e-03
## 617 20:17915546:G:A 8.913829e-03
## 392 16:4021734:A:G 8.923359e-03
## 711 3:55411763:A:G 8.950495e-03
## 697 3:187438522:C:T 8.955303e-03
## 760 4:156620217:G:A 8.971992e-03
## 1042 8:76230301:G:T 8.997585e-03
## 310 14:24935420:C:T 9.022841e-03
## 1071 9:137520289:C:A 9.027615e-03
## 72 1:26803430:T:C 9.050501e-03
## 924 6:44539761:G:A 9.062158e-03
## 161 11:122845075:A:C 9.066913e-03
## 100 1:95787223:T:C 9.135859e-03
## 749 4:130128403:C:T 9.144696e-03
## 496 19:30294991:G:T 9.159304e-03
## 150 10:8701219:G:A 9.162155e-03
## 198 11:69379161:C:A 9.188293e-03
## 793 4:92187079:A:G 9.191245e-03
## 592 2:37517566:A:G 9.191601e-03
## 979 7:28205303:T:C 9.193631e-03
## 233 12:31115718:C:T 9.194909e-03
## 367 15:71772995:C:T 9.195699e-03
## 89 1:67934430:G:A 9.203986e-03
## 424 17:38596425:G:A 9.217332e-03
## 486 19:11206575:G:A 9.245448e-03
## 111 10:105644736:T:G 9.254087e-03
## 825 5:171000733:G:A 9.264757e-03
## 331 14:87626747:C:T 9.286716e-03
## 92 1:81334027:G:A 9.289009e-03
## 306 13:96018966:T:C 9.295669e-03
## 1032 8:3919130:G:A 9.302726e-03
## 573 2:219949184:C:T 9.320609e-03
## 4 1:109817590:G:T 9.322999e-03
## 251 12:65677086:T:C 9.324249e-03
## 696 3:186237034:T:C 9.333554e-03
## 645 20:6699595:T:G 9.366891e-03
## 1043 8:76419046:C:A 9.369349e-03
## 554 2:198688084:A:C 9.384552e-03
## 987 7:91630620:G:T 9.412806e-03
## 823 5:168250903:G:A 9.430376e-03
## 314 14:37144516:C:T 9.432031e-03
## 126 10:22032942:A:G 9.445340e-03
## 29 1:162162334:T:C 9.462020e-03
## 630 20:44608901:G:A 9.465076e-03
## 866 6:109723939:G:A 9.466222e-03
## 471 18:46270114:A:G 9.467056e-03
## 6 1:114448389:C:T 9.468467e-03
## 454 17:74256913:C:T 9.469222e-03
## 481 18:76534379:A:G 9.477307e-03
## 928 6:48279464:C:T 9.477863e-03
## 905 6:2167225:T:C 9.481914e-03
## 664 3:114214611:C:A 9.498315e-03
## 709 3:49913705:T:C 9.504319e-03
## 1053 9:109181911:C:T 9.508507e-03
## 199 11:69933717:T:C 9.517334e-03
## 742 4:111894848:T:G 9.520491e-03
## 882 6:144079629:T:G 9.559551e-03
## 927 6:47487762:A:G 9.575278e-03
## 81 1:50582172:A:C 9.583545e-03
## 508 19:48188809:G:T 9.630748e-03
## 60 1:227750068:G:A 9.743950e-03
## 184 11:61552680:G:T 9.749691e-03
## 479 18:74325253:G:T 9.806590e-03
## 419 17:35469925:T:C 9.881380e-03
## 937 6:76164589:C:A 9.908031e-03
## 843 5:39426020:T:C 9.921533e-03
## 313 14:37128564:C:A 9.960913e-03
## 721 3:68622366:G:A 9.976537e-03
## 355 15:60591309:G:A 1.002360e-02
## 359 15:63039118:C:T 1.003094e-02
## 166 11:13277961:C:T 1.004893e-02
## 1040 8:72851034:C:T 1.004989e-02
## 255 12:69827658:G:T 1.006826e-02
## 455 17:7440584:T:C 1.007167e-02
## 741 4:111856631:T:G 1.007463e-02
## 54 1:222045446:G:T 1.008380e-02
## 493 19:18570035:G:A 1.010112e-02
## 726 4:103687208:A:G 1.012240e-02
## 956 7:124396645:G:A 1.012968e-02
## 737 4:111732536:C:T 1.015355e-02
## 256 12:69936101:G:A 1.015729e-02
## 712 3:56580015:A:G 1.021257e-02
## 66 1:243618317:G:T 1.024158e-02
## 282 13:40817751:T:C 1.027394e-02
## 149 10:82251514:T:C 1.029654e-02
## 665 3:114679536:G:A 1.031544e-02
## 925 6:45244415:T:C 1.032591e-02
## 25 1:154991389:T:C 1.033166e-02
## 1047 9:100931113:C:T 1.034209e-02
## 296 13:75060578:C:T 1.034303e-02
## 976 7:26605282:C:T 1.036438e-02
## 33 1:172083881:G:A 1.036669e-02
## 534 2:169707428:C:T 1.037591e-02
## 441 17:62355081:T:C 1.038435e-02
## 357 15:61408362:G:A 1.040544e-02
## 769 4:21298221:G:T 1.042375e-02
## 18 1:151551320:G:A 1.046542e-02
## 938 6:76173832:C:T 1.047086e-02
## 91 1:78623626:C:T 1.049803e-02
## 446 17:68075324:T:C 1.052487e-02
## 143 10:76861680:T:C 1.055810e-02
## 235 12:4086785:A:C 1.056178e-02
## 385 15:99633448:G:A 1.057109e-02
## 325 14:69034682:C:T 1.060634e-02
## 623 20:35551198:A:C 1.063983e-02
## 747 4:120901336:A:G 1.066956e-02
## 991 7:96013072:A:G 1.068594e-02
## 110 10:105577409:G:A 1.070371e-02
## 863 5:92989489:A:C 1.073510e-02
## 955 7:120777619:G:A 1.073962e-02
## 121 10:123396806:G:T 1.075008e-02
## 113 10:114726843:G:A 1.076247e-02
## 715 3:60348977:A:G 1.078855e-02
## 758 4:156408950:G:A 1.079500e-02
## 265 12:97439396:A:G 1.079903e-02
## 465 18:28750750:A:G 1.080492e-02
## 124 10:127673877:C:T 1.082096e-02
## 422 17:36929578:C:T 1.083062e-02
## 24 1:154911689:G:A 1.083147e-02
## 1007 8:128407443:A:G 1.083156e-02
## 838 5:32771938:A:G 1.083702e-02
## 401 16:73068515:T:C 1.085898e-02
## 906 6:25520094:G:A 1.089152e-02
## 335 14:92926952:G:T 1.090228e-02
## 449 17:70018911:T:G 1.091823e-02
## 1001 8:124020903:C:T 1.093879e-02
## 977 7:2795957:A:G 1.094235e-02
## 656 21:36119111:T:G 1.095080e-02
## 632 20:48634821:G:A 1.095898e-02
## 46 1:208024820:C:T 1.097601e-02
## 127 10:26043071:A:G 1.098880e-02
## 48 1:214627419:C:T 1.099090e-02
## 157 11:114337188:C:T 1.099586e-02
## 165 11:12924265:G:A 1.101066e-02
## 370 15:73650824:G:A 1.102925e-02
## 193 11:67090571:C:T 1.102946e-02
## 1008 8:128413305:G:T 1.109013e-02
## 1016 8:141060747:G:A 1.111591e-02
## 120 10:123334930:G:A 1.116426e-02
## 476 18:55952139:T:C 1.119704e-02
## 172 11:2171601:T:C 1.120270e-02
## 564 2:216291359:C:T 1.122330e-02
## 298 13:76374795:A:G 1.123226e-02
## 7 1:118492052:T:C 1.127590e-02
## 603 2:48642541:G:A 1.128033e-02
## 67 1:247696047:G:A 1.128737e-02
## 985 7:73304636:C:T 1.130562e-02
## 800 5:11114390:G:A 1.130744e-02
## 1039 8:57155598:G:T 1.131450e-02
## 42 1:196316761:C:T 1.133724e-02
## 26 1:159892088:G:A 1.135498e-02
## 1044 8:78153243:C:T 1.142801e-02
## 1036 8:5519379:C:T 1.143523e-02
## 824 5:170875097:A:C 1.143862e-02
## 291 13:59938758:C:A 1.147094e-02
## 615 2:88924622:T:C 1.149292e-02
## 970 7:19248278:G:T 1.153476e-02
## 1050 9:107706188:C:T 1.154913e-02
## 297 13:75624993:T:C 1.157532e-02
## 718 3:65729827:T:C 1.164382e-02
## 146 10:80930698:A:G 1.165836e-02
## 194 11:68918985:A:C 1.166334e-02
## 643 20:6469596:G:A 1.171731e-02
## 1006 8:128352234:G:A 1.172391e-02
## 1015 8:135650483:G:A 1.173783e-02
## 931 6:51834020:C:A 1.179562e-02
## 724 3:9084534:C:T 1.179996e-02
## 641 20:62269750:C:T 1.180477e-02
## 815 5:14568126:G:A 1.183752e-02
## 189 11:65336819:C:T 1.184196e-02
## 5 1:113202571:G:A 1.188505e-02
## 981 7:39418734:C:T 1.190380e-02
## 223 12:13867668:A:G 1.194540e-02
## 659 21:39966133:A:G 1.194832e-02
## 139 10:5898224:G:A 1.195031e-02
## 954 7:116218086:C:T 1.196339e-02
## 826 5:171203438:C:A 1.197505e-02
## 450 17:73368985:A:G 1.200339e-02
## 719 3:66835382:A:C 1.206122e-02
## 229 12:20857467:A:C 1.207828e-02
## 228 12:20536371:G:A 1.212957e-02
## 958 7:131059056:G:A 1.213246e-02
## 485 19:11188247:G:A 1.220366e-02
## 740 4:111832768:C:T 1.220980e-02
## 209 12:107338631:C:T 1.225495e-02
## 507 19:46190268:G:A 1.227147e-02
## 339 15:100246718:C:T 1.228084e-02
## 757 4:152178631:T:C 1.230382e-02
## 74 1:37962756:A:C 1.230659e-02
## 469 18:42399590:A:G 1.231503e-02
## 959 7:135045786:T:G 1.234218e-02
## 459 17:77781725:A:G 1.234613e-02
## 540 2:177598083:A:G 1.235597e-02
## 1034 8:49409929:T:C 1.238170e-02
## 56 1:22270907:C:T 1.238729e-02
## 792 4:88194170:A:G 1.240545e-02
## 478 18:57838401:G:A 1.243658e-02
## 585 2:25483121:G:A 1.244669e-02
## 896 6:161183526:G:A 1.246523e-02
## 751 4:139284343:G:A 1.248735e-02
## 517 2:109047190:A:G 1.249121e-02
## 326 14:70557927:T:G 1.253139e-02
## 203 11:84308799:A:G 1.253455e-02
## 903 6:19841493:T:C 1.255905e-02
## 427 17:42060631:A:G 1.258411e-02
## 988 7:92244422:C:T 1.265216e-02
## 635 20:54842378:A:G 1.271940e-02
## 619 20:21180259:C:T 1.275604e-02
## 207 12:104290997:G:A 1.277050e-02
## 532 2:1645673:C:T 1.282499e-02
## 948 6:97060124:G:A 1.283094e-02
## 304 13:92962685:C:A 1.286869e-02
## 418 17:30239698:T:C 1.289173e-02
## 893 6:161010118:A:G 1.289593e-02
## 188 11:6221249:A:G 1.290349e-02
## 601 2:44702978:T:C 1.292866e-02
## 535 2:172152646:C:T 1.292972e-02
## 1025 8:24092500:C:T 1.293100e-02
## 226 12:19142169:C:T 1.293683e-02
## 490 19:15982760:T:C 1.295277e-02
## 404 16:80650805:A:G 1.296747e-02
## 951 7:116058859:A:G 1.297931e-02
## 684 3:172163449:G:A 1.300390e-02
## 916 6:35614026:T:C 1.301458e-02
## 874 6:122452329:T:C 1.303528e-02
## 399 16:68820946:G:A 1.304216e-02
## 287 13:51084173:T:G 1.305639e-02
## 971 7:19642100:G:T 1.306853e-02
## 821 5:165716407:A:G 1.310571e-02
## 495 19:2176403:G:A 1.311035e-02
## 387 16:1904506:G:A 1.312013e-02
## 1086 9:94228244:G:A 1.312765e-02
## 847 5:56053535:T:C 1.314245e-02
## 962 7:143110762:G:A 1.315933e-02
## 373 15:78691740:G:A 1.319388e-02
## 586 2:27741237:T:C 1.319797e-02
## 530 2:151374518:G:T 1.319927e-02
## 853 5:73964660:A:G 1.322136e-02
## 185 11:61565908:T:C 1.322912e-02
## 259 12:90008959:A:G 1.323773e-02
## 136 10:4965434:A:G 1.332924e-02
## 682 3:169492101:C:T 1.333032e-02
## 278 13:30172751:G:T 1.335316e-02
## 365 15:70048157:A:C 1.336569e-02
## 59 1:227254380:C:T 1.337241e-02
## 1074 9:16787670:C:T 1.337328e-02
## 431 17:47060322:C:A 1.339406e-02
## 168 11:18218678:G:A 1.339543e-02
## 162 11:125497247:T:C 1.342301e-02
## 472 18:4631497:G:A 1.345872e-02
## 897 6:162110497:G:A 1.349473e-02
## 702 3:3647609:T:C 1.356001e-02
## 1004 8:128095156:A:G 1.358868e-02
## 583 2:25139318:C:T 1.360712e-02
## 694 3:185623175:C:T 1.364671e-02
## 828 5:172669771:A:G 1.364780e-02
## 133 10:4011350:A:G 1.366490e-02
## 818 5:157820602:C:T 1.367551e-02
## 714 3:58013573:C:T 1.368154e-02
## 51 1:217718789:C:T 1.372567e-02
## 746 4:120723085:C:T 1.373840e-02
## 461 18:13094132:G:A 1.376425e-02
## 579 2:234068476:C:T 1.379378e-02
## 589 2:33405151:T:C 1.379976e-02
## 423 17:38155798:G:A 1.383131e-02
## 625 20:39780932:G:A 1.387374e-02
## 851 5:65942372:A:G 1.389160e-02
## 855 5:77505876:C:T 1.389343e-02
## 451 17:7371932:C:T 1.389404e-02
## 57 1:224630695:C:T 1.390536e-02
## 55 1:222164948:A:G 1.392256e-02
## 93 1:84172192:A:G 1.393689e-02
## 599 2:43629612:A:G 1.401268e-02
## 484 19:10991827:G:A 1.403336e-02
## 550 2:191604422:G:A 1.404597e-02
## 109 10:105340059:A:G 1.404734e-02
## 68 1:25032168:G:A 1.404975e-02
## 996 8:116637685:T:C 1.406226e-02
## 489 19:15655168:C:T 1.407421e-02
## 596 2:43167878:A:C 1.408267e-02
## 47 1:212237798:G:A 1.409412e-02
## 37 1:183081194:A:C 1.410195e-02
## 227 12:19531205:T:C 1.412186e-02
## 934 6:72203005:G:A 1.412517e-02
## 403 16:80025366:G:A 1.414803e-02
## 377 15:89356832:G:T 1.415665e-02
## 725 3:98715823:G:A 1.416039e-02
## 420 17:36074979:G:A 1.416924e-02
## 409 16:86688976:C:T 1.417977e-02
## 397 16:5669598:G:T 1.421911e-02
## 784 4:63248138:A:G 1.422461e-02
## 498 19:33532300:C:T 1.423331e-02
## 960 7:139463565:A:G 1.426672e-02
## 237 12:42796300:A:G 1.427130e-02
## 887 6:152345162:G:A 1.429575e-02
## 70 1:25753638:C:T 1.432508e-02
## 246 12:576984:C:T 1.442672e-02
## 393 16:50259483:A:G 1.447311e-02
## 69 1:25044111:C:T 1.447776e-02
## 52 1:218027794:T:C 1.447915e-02
## 689 3:18447100:A:G 1.452958e-02
## 787 4:75515054:A:G 1.453922e-02
## 713 3:56585975:T:C 1.460277e-02
## 672 3:138119952:T:C 1.461637e-02
## 112 10:105671527:T:C 1.461992e-02
## 658 21:37395301:T:C 1.462337e-02
## 819 5:158220193:A:G 1.464124e-02
## 794 5:107343923:T:C 1.471446e-02
## 346 15:36039864:T:C 1.471914e-02
## 381 15:94028149:T:C 1.472802e-02
## 78 1:41560235:C:T 1.474985e-02
## 648 21:19040005:T:C 1.475061e-02
## 761 4:16035795:A:C 1.477155e-02
## 219 12:122494809:T:C 1.479527e-02
## 519 2:127423075:G:A 1.481970e-02
## 842 5:38768482:A:G 1.485444e-02
## 894 6:161097871:A:G 1.490863e-02
## 213 12:115094260:A:G 1.491533e-02
## 1051 9:108304500:C:A 1.493027e-02
## 700 3:191441507:C:T 1.495665e-02
## 649 21:24339568:C:T 1.496849e-02
## 321 14:61008596:T:C 1.496937e-02
## 912 6:33628863:A:G 1.498044e-02
## 1019 8:143427682:G:A 1.500710e-02
## 773 4:38688362:C:T 1.501909e-02
## 1030 8:29509616:A:C 1.502187e-02
## 1013 8:135612595:A:G 1.506108e-02
## 406 16:82184201:A:C 1.506595e-02
## 556 2:198952637:C:T 1.510185e-02
## 204 11:84524271:C:A 1.512725e-02
## 336 14:93104072:T:C 1.516039e-02
## 856 5:77837789:G:A 1.517411e-02
## 743 4:114741562:A:G 1.518717e-02
## 904 6:21163919:T:C 1.521166e-02
## 545 2:18574952:T:C 1.522431e-02
## 258 12:87301493:A:G 1.524856e-02
## 173 11:2286243:G:A 1.526838e-02
## 217 12:11855773:G:A 1.527144e-02
## 104 10:102684380:C:A 1.530936e-02
## 407 16:84987679:C:A 1.536023e-02
## 1020 8:16816092:C:T 1.536756e-02
## 232 12:29496991:C:T 1.538022e-02
## 876 6:122708686:G:A 1.538549e-02
## 745 4:120509972:G:A 1.539944e-02
## 122 10:12488714:T:G 1.544903e-02
## 642 20:6404281:C:A 1.545422e-02
## 305 13:94738136:A:G 1.553537e-02
## 518 2:124482626:C:T 1.556098e-02
## 546 2:190387487:C:T 1.557972e-02
## 211 12:114797093:G:T 1.559191e-02
## 836 5:32567732:C:T 1.559321e-02
## 1087 9:94258836:A:C 1.559834e-02
## 932 6:70987200:G:T 1.560730e-02
## 430 17:46974734:T:C 1.562278e-02
## 972 7:20292134:G:A 1.571971e-02
## 850 5:64674446:C:T 1.574064e-02
## 830 5:173377817:G:A 1.576692e-02
## 614 2:86390714:G:A 1.577486e-02
## 1023 8:17893093:A:C 1.578647e-02
## 1011 8:130384476:T:C 1.583051e-02
## 474 18:48142854:T:G 1.583490e-02
## 344 15:32994756:T:C 1.583591e-02
## 201 11:75276178:A:C 1.583740e-02
## 944 6:82128386:T:C 1.584208e-02
## 947 6:85448103:C:T 1.584578e-02
## 533 2:165558252:T:C 1.584590e-02
## 14 1:145644984:C:T 1.586427e-02
## 651 21:28704127:T:C 1.587112e-02
## 415 17:2132324:T:C 1.597010e-02
## 435 17:54839652:A:G 1.598190e-02
## 11 1:120303055:C:T 1.605306e-02
## 949 7:100523:T:C 1.606613e-02
## 610 2:70172587:G:A 1.607137e-02
## 565 2:217920769:G:T 1.611057e-02
## 872 6:118575602:G:A 1.620775e-02
## 921 6:41536427:C:T 1.621498e-02
## 1088 9:95387983:C:T 1.623925e-02
## 588 2:33361425:A:G 1.624994e-02
## 884 6:149755695:A:G 1.629217e-02
## 857 5:80006558:C:T 1.631708e-02
## 920 6:41290458:G:A 1.637196e-02
## 356 15:61294826:T:C 1.643085e-02
## 1070 9:136149830:G:A 1.645811e-02
## 552 2:192587204:T:C 1.648545e-02
## 240 12:50901882:T:C 1.649446e-02
## 391 16:30927509:C:T 1.650447e-02
## 1091 9:97571249:T:C 1.650950e-02
## 269 13:111114176:C:A 1.652651e-02
## 262 12:94126925:T:G 1.653120e-02
## 1062 9:118305735:A:G 1.658174e-02
## 738 4:111745599:A:G 1.663135e-02
## 116 10:119621825:A:G 1.664866e-02
## 890 6:157153670:G:A 1.665465e-02
## 366 15:70364352:G:A 1.666770e-02
## 192 11:66826160:C:T 1.680833e-02
## 1046 9:100482976:C:A 1.702452e-02
## 206 12:103062597:A:G 1.705548e-02
## 86 1:64098065:T:C 1.705908e-02
## 71 1:26450009:C:T 1.706257e-02
## 634 20:52156847:T:G 1.709616e-02
## 854 5:75038431:T:G 1.711731e-02
## 622 20:34434916:G:A 1.714055e-02
## 986 7:77486676:T:C 1.716514e-02
## 408 16:86417890:C:T 1.719784e-02
## 1063 9:118640119:C:T 1.719860e-02
## 940 6:7792947:T:C 1.723114e-02
## 671 3:137461111:A:G 1.727392e-02
## 772 4:37926978:C:T 1.729296e-02
## 662 3:106779118:A:G 1.734870e-02
## 786 4:73476014:T:G 1.736845e-02
## 526 2:145231349:A:G 1.739328e-02
## 22 1:154814353:C:T 1.739738e-02
## 97 1:93343282:A:G 1.746693e-02
## 1068 9:130932776:C:T 1.750069e-02
## 766 4:18010384:A:G 1.754151e-02
## 570 2:218688596:T:G 1.759716e-02
## 85 1:59653742:A:G 1.770383e-02
## 354 15:58650355:C:T 1.770944e-02
## 990 7:92659160:C:T 1.771433e-02
## 978 7:27976563:G:A 1.775104e-02
## 568 2:218296508:C:T 1.777715e-02
## 964 7:145720467:A:G 1.778719e-02
## 167 11:14404825:C:T 1.786384e-02
## 358 15:62379971:A:G 1.788374e-02
## 410 16:89875710:A:C 1.791630e-02
## 748 4:12963574:G:A 1.793358e-02
## 31 1:170346164:C:T 1.793654e-02
## 717 3:63967900:A:G 1.794274e-02
## 178 11:35972672:G:A 1.807000e-02
## 834 5:179730843:A:G 1.809441e-02
## 963 7:144074929:G:A 1.809577e-02
## 881 6:142703877:G:A 1.810963e-02
## 215 12:116391296:T:C 1.811836e-02
## 353 15:54455757:T:C 1.812501e-02
## 1012 8:130723728:A:G 1.816356e-02
## 123 10:126383088:C:A 1.820138e-02
## 106 10:104604916:T:G 1.822156e-02
## 505 19:45415713:G:A 1.827054e-02
## 1083 9:86661654:C:T 1.827304e-02
## 9 1:118854698:G:A 1.829201e-02
## 946 6:84277037:C:A 1.829773e-02
## 626 20:40534202:G:A 1.833764e-02
## 129 10:27890831:A:G 1.839150e-02
## 553 2:198077877:A:G 1.843233e-02
## 512 19:7244233:A:C 1.844957e-02
## 1009 8:128424792:A:G 1.851813e-02
## 595 2:42462930:A:G 1.855264e-02
## 372 15:78686823:A:G 1.866659e-02
## 2 1:10566215:A:G 1.871633e-02
## 473 18:46657358:C:T 1.872113e-02
## 340 15:101632867:G:A 1.877734e-02
## 230 12:24784139:G:A 1.885831e-02
## 361 15:66867861:A:G 1.889948e-02
## 597 2:43463637:T:C 1.893412e-02
## 275 13:23371474:T:C 1.898917e-02
## 807 5:132949899:A:G 1.899640e-02
## 443 17:6474096:C:T 1.899675e-02
## 53 1:218975475:G:A 1.900778e-02
## 247 12:58256714:T:G 1.905047e-02
## 432 17:47402807:C:T 1.906915e-02
## 609 2:68495002:G:A 1.907063e-02
## 330 14:75038472:T:C 1.907516e-02
## 205 11:85867875:A:G 1.910143e-02
## 65 1:243471192:A:G 1.913178e-02
## 1045 8:93643870:C:T 1.914344e-02
## 183 11:59923508:A:G 1.914429e-02
## 652 21:30567941:G:T 1.917545e-02
## 284 13:49201040:T:G 1.922188e-02
## 1090 9:97429160:A:G 1.923755e-02
## 1079 9:33026572:T:C 1.930564e-02
## 538 2:174815898:A:C 1.931885e-02
## 762 4:163337191:C:T 1.933543e-02
## 1072 9:138698196:G:A 1.940762e-02
## 704 3:38624343:A:G 1.942865e-02
## 389 16:2248888:T:C 1.943051e-02
## 862 5:90732225:C:T 1.950967e-02
## 868 6:117210052:T:C 1.951065e-02
## 274 13:23329181:G:A 1.954871e-02
## 1052 9:108901049:C:T 1.955231e-02
## 115 10:116828824:T:G 1.963734e-02
## 294 13:66712732:G:A 1.970178e-02
## 1024 8:19493786:T:C 1.970229e-02
## 360 15:63439186:C:T 1.970537e-02
## 870 6:117868051:T:C 1.978912e-02
## 266 13:111015333:G:T 1.983532e-02
## 23 1:154834183:G:A 1.993863e-02
## 270 13:111380701:T:G 2.008756e-02
## 1021 8:17292466:T:G 2.010127e-02
## 293 13:63436800:A:G 2.012624e-02
## 243 12:54177548:A:G 2.014170e-02
## 957 7:130753991:C:T 2.015067e-02
## 705 3:38795555:A:G 2.016510e-02
## 594 2:38276549:A:G 2.017851e-02
## 933 6:7213016:G:A 2.026708e-02
## 62 1:23504795:A:G 2.031655e-02
## 989 7:92300587:G:A 2.033230e-02
## 3 1:108018847:C:T 2.040417e-02
## 177 11:30347927:A:G 2.040868e-02
## 654 21:33325270:A:G 2.043783e-02
## 604 2:56089540:A:G 2.044150e-02
## 969 7:19049388:G:A 2.044262e-02
## 891 6:158722034:A:G 2.049026e-02
## 861 5:90151589:G:A 2.049371e-02
## 673 3:141105570:A:G 2.049497e-02
## 380 15:91512067:G:A 2.050657e-02
## 953 7:116191301:C:A 2.053022e-02
## 613 2:85491365:T:C 2.056072e-02
## 128 10:26512375:C:T 2.059316e-02
## 386 16:14388626:C:T 2.061232e-02
## 491 19:17522309:T:C 2.062193e-02
## 405 16:81574197:T:C 2.064536e-02
## 605 2:56096892:A:G 2.072007e-02
## 1027 8:26855424:C:A 2.072776e-02
## 1048 9:101743336:T:C 2.075985e-02
## 753 4:145574844:A:G 2.077589e-02
## 563 2:216245392:C:T 2.098995e-02
## 436 17:55959944:A:G 2.103415e-02
## 283 13:47112120:C:A 2.105276e-02
## 402 16:75328308:T:C 2.115352e-02
## 1075 9:21801530:A:G 2.132602e-02
## 877 6:126717064:A:G 2.136116e-02
## 612 2:71627539:C:T 2.136588e-02
## 692 3:185488882:G:A 2.144078e-02
## 628 20:4098567:A:G 2.154633e-02
## 1033 8:42324765:T:C 2.165147e-02
## 345 15:33769513:C:T 2.170626e-02
## 371 15:73664310:T:G 2.176644e-02
## 1078 9:31489011:T:C 2.177560e-02
## 898 6:16418562:C:T 2.180897e-02
## 678 3:157806960:T:C 2.183216e-02
## 816 5:153970243:G:A 2.187997e-02
## 39 1:193424164:T:C 2.188539e-02
## 234 12:32972940:C:T 2.191417e-02
## 980 7:32935224:G:T 2.200121e-02
## 559 2:202799924:T:C 2.201945e-02
## 378 15:89395626:A:G 2.203182e-02
## 529 2:145801113:G:A 2.205118e-02
## 926 6:47475022:T:G 2.206290e-02
## 463 18:24570667:T:G 2.206322e-02
## 318 14:54560018:T:C 2.210907e-02
## 660 21:40426768:A:G 2.211514e-02
## 1094 9:98318926:C:T 2.211695e-02
## 782 4:57838583:G:T 2.212397e-02
## 537 2:174212894:G:A 2.213075e-02
## 1029 8:27527995:C:T 2.213876e-02
## 637 20:57472174:C:T 2.213893e-02
## 820 5:158244083:C:T 2.214279e-02
## 1057 9:110895353:C:T 2.215313e-02
## 581 2:24247514:A:G 2.215940e-02
## 179 11:45274894:C:T 2.216225e-02
## 813 5:142792484:G:T 2.217652e-02
## 155 11:10657739:C:T 2.223775e-02
## 560 2:203194256:A:C 2.227200e-02
## 655 21:35659776:G:A 2.227678e-02
## 195 11:68994667:A:G 2.228755e-02
## 322 14:63928710:T:C 2.232707e-02
## 252 12:66359752:C:A 2.238603e-02
## 506 19:45422846:G:A 2.243868e-02
## 375 15:81013037:G:A 2.244044e-02
## 846 5:53361606:T:C 2.247655e-02
## 516 2:108643385:A:G 2.258235e-02
## 190 11:65390803:G:A 2.269053e-02
## 163 11:12679245:T:G 2.271505e-02
## 202 11:80644946:T:G 2.273537e-02
## 497 19:31047269:A:G 2.283462e-02
## 1041 8:76040583:C:T 2.293474e-02
## 878 6:12903957:A:G 2.297638e-02
## 501 19:41922352:A:G 2.299547e-02
## 900 6:168834623:C:T 2.301010e-02
## 80 1:49070646:T:C 2.310176e-02
## 492 19:17844894:A:G 2.317985e-02
## 263 12:94880742:G:A 2.318080e-02
## 804 5:127868199:G:A 2.333346e-02
## 369 15:73370467:G:A 2.337196e-02
## 447 17:69108753:G:T 2.340381e-02
## 860 5:88327782:T:G 2.341532e-02
## 771 4:36327093:G:A 2.341975e-02
## 611 2:71584485:G:A 2.342075e-02
## 910 6:32373185:T:C 2.345255e-02
## 549 2:191437024:A:G 2.354694e-02
## 582 2:25040082:C:T 2.365119e-02
## 27 1:160318940:A:G 2.367938e-02
## 973 7:20381674:G:A 2.370621e-02
## 809 5:135678996:T:C 2.372241e-02
## 1 1:103519589:G:A 2.387109e-02
## 914 6:34773727:C:T 2.390860e-02
## 320 14:59688820:G:A 2.392682e-02
## 966 7:150508720:C:T 2.398928e-02
## 1077 9:27982137:A:G 2.412054e-02
## 140 10:70196580:T:C 2.421556e-02
## 28 1:160399586:T:C 2.432670e-02
## 144 10:80819132:A:G 2.432806e-02
## 316 14:54387873:C:T 2.435820e-02
## 445 17:66303352:G:A 2.443567e-02
## 685 3:173179298:C:T 2.458142e-02
## 452 17:73908566:A:C 2.463523e-02
## 1092 9:97713459:C:A 2.473967e-02
## 12 1:13885203:C:T 2.481426e-02
## 272 13:22319948:T:C 2.481499e-02
## 458 17:76718842:C:T 2.498174e-02
## 504 19:44286513:A:G 2.501058e-02
## 499 19:38346685:A:G 2.501788e-02
## 285 13:50469913:C:T 2.505273e-02
## 998 8:117563532:G:A 2.514788e-02
## 302 13:81550449:T:C 2.521085e-02
## 584 2:25463483:G:A 2.525806e-02
## 551 2:191832662:T:C 2.538334e-02
## 679 3:158285175:T:C 2.544736e-02
## 608 2:57218797:G:T 2.552728e-02
## 148 10:82170856:A:G 2.555518e-02
## 453 17:73949045:C:T 2.556254e-02
## 1014 8:135612745:A:G 2.557238e-02
## 19 1:151763246:G:A 2.566957e-02
## 277 13:29022645:C:T 2.567455e-02
## 983 7:46201355:A:G 2.571275e-02
## 398 16:57548744:T:C 2.584089e-02
## 520 2:128050134:T:C 2.586018e-02
## 118 10:121129797:G:A 2.586723e-02
## 32 1:170569689:G:A 2.591555e-02
## 156 11:111171709:C:A 2.595510e-02
## 250 12:64971147:G:A 2.613864e-02
## 796 5:107975671:A:G 2.618895e-02
## 621 20:34025756:A:G 2.620091e-02
## 169 11:1909006:T:C 2.623693e-02
## 693 3:185548683:G:A 2.628930e-02
## 952 7:116154015:T:G 2.630024e-02
## 710 3:53134098:T:C 2.637836e-02
## 620 20:33718706:G:A 2.667249e-02
## 997 8:117209548:A:G 2.667751e-02
## 117 10:120738290:C:T 2.670669e-02
## 280 13:33143406:G:A 2.672287e-02
## 646 20:8626271:C:A 2.684549e-02
## 727 4:103936001:G:A 2.700609e-02
## 731 4:107515773:A:G 2.710300e-02
## 442 17:64473021:A:G 2.726536e-02
## 859 5:81996389:A:G 2.735094e-02
## 488 19:1384690:G:A 2.736502e-02
## 683 3:172046933:A:G 2.737834e-02
## 789 4:82155568:C:A 2.738706e-02
## 915 6:35246903:G:A 2.745500e-02
## 279 13:33057345:G:T 2.755676e-02
## 241 12:51155663:C:T 2.762854e-02
## 732 4:1078124:G:A 2.769558e-02
## 131 10:32082658:T:C 2.788147e-02
## 35 1:176794066:G:A 2.795765e-02
## 974 7:23475919:G:A 2.802421e-02
## 317 14:54410919:T:C 2.806008e-02
## 767 4:181909685:G:A 2.832432e-02
## 1061 9:118169080:A:G 2.833457e-02
## 238 12:4368352:T:C 2.839215e-02
## 396 16:53842908:G:A 2.840857e-02
## 624 20:38552078:T:C 2.864745e-02
## 1018 8:141992778:T:C 2.865140e-02
## 1017 8:141729463:C:T 2.881498e-02
## 879 6:134173151:A:G 2.887869e-02
## 412 17:1054071:T:C 2.899031e-02
## 75 1:38386727:G:A 2.899863e-02
## 795 5:107566842:C:T 2.911059e-02
## 102 10:101805442:C:T 2.926966e-02
## 159 11:117283676:T:C 2.937196e-02
## 323 14:65567515:G:T 2.937272e-02
## 869 6:117490664:T:C 2.937322e-02
## 566 2:217935116:T:C 2.953341e-02
## 108 10:105324774:C:T 2.955258e-02
## 295 13:73957681:A:G 2.955401e-02
## 528 2:145645775:G:A 2.980645e-02
## 841 5:33989518:C:T 2.987911e-02
## 752 4:145565826:C:T 3.046719e-02
## 640 20:60958269:T:C 3.073592e-02
## 363 15:67446831:G:A 3.083248e-02
## 16 1:149938898:T:C 3.104603e-02
## 852 5:67599656:T:G 3.122608e-02
## 776 4:41261500:T:C 3.129284e-02
## 421 17:36101156:T:C 3.134972e-02
## 577 2:232322779:G:A 3.146841e-02
## 448 17:69923355:C:T 3.150043e-02
## 303 13:87001719:T:C 3.152685e-02
## 587 2:28635740:T:C 3.163465e-02
## 218 12:121416988:A:G 3.163669e-02
## 522 2:136187345:C:T 3.169046e-02
## 888 6:152437016:T:C 3.183736e-02
## 908 6:28747902:T:C 3.196243e-02
## 130 10:30323892:G:A 3.200226e-02
## 840 5:33499437:A:G 3.207782e-02
## 961 7:139942304:G:A 3.210554e-02
## 239 12:4942953:T:G 3.219421e-02
## 524 2:142746083:G:A 3.231244e-02
## 708 3:4742276:A:G 3.249171e-02
## 873 6:121781390:T:C 3.263183e-02
## 1031 8:36858483:A:G 3.279838e-02
## 668 3:134379752:C:T 3.281386e-02
## 909 6:31592820:C:T 3.294651e-02
## 145 10:80845746:G:A 3.320756e-02
## 231 12:26228883:G:A 3.352449e-02
## 509 19:5146582:A:G 3.355730e-02
## 1049 9:105754656:T:C 3.364954e-02
## 191 11:65583066:G:T 3.380290e-02
## 10 1:119491784:A:G 3.438420e-02
## 483 19:10803395:A:G 3.464492e-02
## 379 15:89851580:G:A 3.475654e-02
## 803 5:127696022:A:C 3.477973e-02
## 562 2:208526140:T:C 3.485917e-02
## 1028 8:27195121:T:C 3.499621e-02
## 76 1:38559352:A:G 3.500686e-02
## 892 6:160833664:C:T 3.506576e-02
## 829 5:172994624:G:A 3.535245e-02
## 319 14:55245149:A:G 3.542141e-02
## 889 6:156587831:T:C 3.559188e-02
## 867 6:116451442:A:G 3.602886e-02
## 222 12:13075379:G:A 3.609164e-02
## 688 3:183355405:A:G 3.614213e-02
## 707 3:46600520:A:G 3.646873e-02
## 729 4:106084778:C:T 3.651783e-02
## 20 1:154404406:A:C 3.670752e-02
## 730 4:106216205:T:C 3.713331e-02
## 706 3:41243742:T:G 3.730801e-02
## 1093 9:98265780:C:T 3.781516e-02
## 382 15:95312071:T:C 3.847165e-02
## 434 17:53209774:A:C 3.849144e-02
## 805 5:1297488:C:T 3.849828e-02
## 513 2:101576199:G:A 3.900736e-02
## 695 3:186004526:T:C 3.933196e-02
## 374 15:79121776:C:T 3.986323e-02
## 844 5:4012694:G:A 3.986729e-02
## 138 10:5886734:C:T 3.987780e-02
## 639 20:60920887:G:A 4.081819e-02
## 885 6:151942194:T:G 4.149024e-02
## 271 13:21570246:T:G 4.210229e-02
## 631 20:48316027:C:T 4.369981e-02
## 907 6:26200677:A:G 4.507931e-02
## 225 12:18040964:A:C 4.572318e-02
## 728 4:106061534:C:A 4.717647e-02
## 734 4:111586468:A:G 4.918761e-02
## 224 12:14520701:A:G 4.978408e-02
## 806 5:131585958:A:G 5.031346e-02
## 514 2:10178479:C:T 5.415964e-02
## 312 14:32981484:G:A 5.612800e-02
## 736 4:111706275:A:G 5.732161e-02
hist(temp.locus_cohort.maxdeltaaf$maxdeltaaaf)
# On the basis of this histogram, set a max delta aaf threshold of
# 4%.
mean(temp.locus_cohort.maxdeltaaf$maxdeltaaaf < 0.04)
## [1] 0.9881279
temp.sel_loci = temp.locus_cohort.maxdeltaaf$vid[temp.locus_cohort.maxdeltaaf$maxdeltaaaf < 0.04]
afs = afs[afs$vid %in% temp.sel_loci,]
models = models[models$vid %in% temp.sel_loci,]
scores = daply(expand.grid(cohort = cohorts.sel, model = models.sel), .(model, cohort), function(cm)
cohort_mean_score(afs.nomissing[afs.nomissing$cohort == cm$cohort,], models[models$id == cm$model,]), .progress = "text")
set.seed(314159)
boot_scores = daply(expand.grid(cohort = cohorts.sel, model = models.sel), .(model, cohort), function(cm)
cohort_mean_score_boot(afs.nomissing[afs.nomissing$cohort == cm$cohort,], models[models$id == cm$model,], B = 5000, drawmissing = TRUE), .progress = "text")
# Calculate all approximate tests, but obviously we will discard most
# of the outputs.
set.seed(314159)
boot_scores_tests = boot_tests(boot_scores, B = 100000)
saveRDS(list(
scores = scores,
boot = boot_scores,
tests = boot_scores_tests), file = "common_variants.rds")
library(gplots)
##
## Attaching package: 'gplots'
## The following object is masked from 'package:stats':
##
## lowess
library(viridis)
## Loading required package: viridisLite
heatmap.2(scores, trace = "none", scale = "row", dendrogram = "none", Colv = FALSE, Rowv = FALSE, col = viridis(100), margin = c(7, 25))
library(ggplot2)
ggplot(adply(boot_scores, c(1,2), function(x) data.frame(boot_score = x)), aes(x = boot_score, col = cohort)) + stat_ecdf() + facet_wrap(~ model, scales = "free") + theme_bw()